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IMPACT OF DEPENDENCE ON SOME MULTIVARIATE RISK INDICATORS 


V. MAUME-DESCHAMPS, D. RULLIERE, AND K. SAID 


Abstract. The minimization of some multivariate risk indicators may be used as an allocation 
method, as proposed in Cenac et al. [3]. The aim of capital allocation is to choose a point in 
a simplex, according to a given criterion. In a previous paper m we proved that the proposed 
allocation technique satisfies a set of coherence axioms. In the present one, we study the properties 
and asymptotic behavior of the allocation for some distribution models. We analyze also the impact 
of the dependence structure on the allocation using some copulas. 


Introduction 

Natural phenomena, financial events and risks are usually modeled through random vectors or 
processes. In these fields, considering dependencies between random variables is necessary. In 
actuarial science, this issue has led to the development of multivariate risk theories. Indeed, an 
insurance company is generally exposed to several risks which cannot be assumed to be indepen¬ 
dent. It is therefore necessary to adopt a multivariate approach that takes into account both the 
marginal structures of risks and their dependence structure. 

Multivariate risk theory is based on dependence modeling and includes multivariate ruin prob¬ 
abilities and multivariate risk measures. In the univariate case, ruin probability has been widely 
studied since the beginning of the 20th century (Lundberg (1903) [12], Cramer (1930) [7|). In mul¬ 
tivariate contexts, several definitions are proposed for ruin probability. Hult and Lindskog (2006) 
[T^ defined a multivariate ruin probability based on a notion of ruin sets. Cai and Li (2007) P| 
defined different finite-time ruin probabilities. Risk measure theory was also enriched by some 
multivariate risk measures definitions. Jouini et al. (2004) [12] defined vector-valued coherent risk 
measures. Cousin and Di Bernardino (2013) [2] introduced some multi-dimensional extensions of 
usual univariate risk measures such as the Value-at-Risk. These measure are multi-dimensional 
valued, and thus cannot take benefit from a full order. Dhaene et al. (2012) introduced of real 
valued family for risk measures for mono-periodic multivariate processes (HU)- In multi-periodic 
context, Cenac et al. (2012) [9] defined new multivariate risk indicators as sum of expected local 
ruin amounts using penalty functions. 

Capital allocation may be seen as a direct application of multivariate risk theory. Allocating 
a capital in actuarial contexts means distributing a positive quantity representing an allocation 
capital in a simplex. For insurance industry the calculation of the regulatory economic capital, 
which is called the Solvency Capital Requirement, is well controlled and its methodology is almost 
imposed by the supervisory authorities of the sector. Nevertheless, the allocation of this capital 
may be considered as an internal exercise for each company, and constitutes a management choice 
whose success is a key factor for firm performance optimization. It can also be seen as an indicator 
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of its good governance. 

Several allocation methods were proposed in the actuarial literature. Some are driven from 
ideas of risk allocation using Euler or Shapley principles (see as examples Tasche (2007) |2I] for 
the Euler allocation principle and Denault(2001) [lU] for the Shapley allocation principle). Others 
are based on the minimization of some ruin probabilities or multivariate risk indicators. In this 
context, Cenac et al. (2012) [H] proposed a capital allocation by minimizing some multivariate 
risk indicators. They presented a numerical optimization algorithm to hnd the optimal allocation 
in some general cases. In [S] properties and asymptotic behavior of the allocation are studied in 
some bivariate cases. The impact of dependence on capital allocation was studied for certain cases 
of Euler’s method by Barges et al. (2009) [2] and Cossette et al. (2012) [3]. 

In this paper, we give a further study on the allocation method by minimizing some multivariate 
risk indicators. The idea of minimizing multivariate risk indicators seems from our point of view 
more suited for the ORS7\0 approach of Solvency 2 European norms. In a recent paper na , we 
show that this method satishes a set of coherence axioms. In the present paper, we study the im¬ 
pact of marginal distributions and dependence structure on the optimal allocation. We generalize 
the results presented in bivariate cases by Cenac et al. (2014) |H]. Moreover, we study the impact 
of the dependence structure on the allocation using some parametric copulas. Our main results 
are explicit formulas for various classes of multivariate processes, as well as the limit behavior of 
the allocation as the capital goes to inhnity. 

The paper is organized as follows. In the hrst section, we recall the allocation method by 
minimizing multivariate risk indicators. The second section is a generalization in higher dimension 
of the results presented for the bivariate case in Cenac et al.(2014) [S]. We present some explicit 
formulas obtained for some particular models, and we discuss the asymptotic behavior of the 
optimal allocation for these models. The impact of the dependence structure on the allocation 
is studied in Section 3 through the analysis of comonotonic cases and some examples of copulas 
(FGM, Marshall-Olkins...). 


1. Optimal allocation 

In this section, we recall one capital allocation principle that consists of minimizing some risk 
indicators. This method was introduced in Cenac et al. (2012) pj. It is based on the minimization 
of some multivariate risk indicators. 

In a multivariate risk framework, we consider an insurance group composed of d branches or 
business lines. We denote by u the initial capital of the group. Let be a vectorial risk process 
Xp = {Xp,... ,Xp), where Xp corresponds to the losses of the business line during the 

p 

period. We denote by the reserve of the line at time p, so: Rp = Uk — '^ X^, where Uk G M"*" 

e=i 

is the initial capital of the business line, then u = Ui + ■ ■ ■ + Ud- 

Cenac et al. (2012) [9] dehned the two following multivariate risk indicators, given penalty 
functions gk- 
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• The indicator /: 


/(mi, . . . , Ud) — X! I X! R^>0} I ’ 

The indicator J: 

J{Ui, . . . , Ud) = X! ( X! 9k{Rp)^{Rk<0}^Ri<0} I ’ 


k=l \p=l 


d / n 


k=l \p=l 


Qk : M“ —>• are convex functions with gk{0) = 0, gk{x) > 0 for a: < 0, k = 1,... ,d. They 

represent the cost that each branch has to pay when it becomes insolvent while the group is solvent 
for the I indicator, or while the group is also insolvent in the case of the J indicator. 

These multivariate risk indicators are a generalization of the mono-periodic indicator family intro¬ 
duced in Dhaene et ah (2012) [TT]. 


The indicator I represents the expected sum of penalty amounts of local ruins, knowing that 
the group remains solvent. In the case of the indicator J, the local ruin severities are taken into 
account only in the case of group insolvency. 


The idea of the optimal allocation is to allocate some capital u by minimizing these indi¬ 
cators. This is hnding an allocation vector {ui,... ,Ud) that minimizes the indicator such as 
u = Ui + ■ ■ ■ + Ud, where u is the initial capital that need to be shared among all branches. In 
our paper im, we explained why this capital allocation method can be considered as economically 
coherent. 


By using optimization stochastic algorithms, one can estimate the minimum of these risk indi¬ 
cators. Cenac et ah (2012) [S] propose a Kiefer-Wolfowitz version of the mirror algorithm as a 
convergent algorithm under general assumptions to hnd optimal allocation minimizing the indictor 
I. This algorithm is effective to solve the optimal allocation problem, especially for a large number 
of business lines, and for allocation over several periods. 

Since the solvency capital requirement (SCR) is calculated for one year time horizon only, it 
seems more practical to focus on the case of mono-periodic allocation. 

The risk Xk corresponds to the losses of the branch during one period. It is a positive ran¬ 
dom variable in our context. We denote by u the initial capital of the hrm, it represents the 
allocation capital and it can be, as an example, the SCR or another investment capital. The 
quantity Uk represents the portion of capital allocated to the branch, then J2i=i Ui = u. The 
simplex = {v = {vi,... ,Vd) E [0, m]'^, = u} is the set of all possible allocations. For 

all i G {!,..., d} we denote a, = then, I]f=i = 1- For {ui, ■ ■ ■ ,Ud) E U^, we dehne also 
the reserve of the business line at the end of the period as = Uk — Xk- The aggregate 
sum of risks is S' = J2f=iXi, and let S“* = for all i E {1,... ,d}. Finally, Fz is the 

distribution function of a random variable Z, Fz is its survival function and fz its density function. 

In order to ensure the existence and uniqueness of the indicators minimum in we suppose 
that for at least one k E {!,... ,d}, gk is strictly convex. We also assume that the joint density 
f{Xk,s) support contains [0,m]^ (see [S]). In this case, the indicators / and J are strictly convex 
and admit a unique minimum. 
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Now we focus on the optimality condition for the indicators / and J. 

For an initial capital u, and an optimal allocation minimizing the multivariate risk indicator I, we 
look for u* G such that: 

I(u*) = inf I(v), V e Ri. 

We assume that the functions gk are differentiable and such that for all k G d}, the 

derivative g'k{uk — X^) admits a moment of order one, and that the random vector (X^, S) has a 
joint density distribution denoted by f(Xk,s)- Under these assumptions, the risk indicators I and 
J are differentiable, and in this case, we can calculate the following gradients: 

^+oo 

(V/(n))i = ^ / gkivk-x)fx^,s{x,u)dx + E[g'i{vi-Xi)l{Xi>vi}Ti-{s<u}] 

^ ^+oo 

and, (VJ(n))i = Xl / 9k{vk - x)fx^:,s{x,u)dx+ E[g'i{vi - Xi)l{Xi>vi}T^{s>u}]- 

k=i 

Using the Lagrange multipliers method, we obtain an optimality condition verihed by the unique 
solution to this optimization problem: 

(1.1) R[g[{ui — Xj)ll|Xi>ni}IL{5<n}] = — Xj)I1{Xj>uj}]1{5<u}] 5 Vj G {1,..., dY- 

A natural choice for penalty functions is the ruin severity: gk{x) = |x|. In that case, and if the joint 
density f{Xk,s) support contains [0 ,m]^, for at least one k E . ,d}, our optimization problem 
has a unique solution. 

In Maume-Deschamps et ah (2015) [I7j, we showed that in the case of penalty functions 
gk{x) = |a;| Vfc G {1,... ,d}, and for continuous random vector (Xi,... ,Xd) such that the joint 
density f(Xk,s) support contains [0 ,m]^, for at least one /c G {1,... ,d}, the optimal allocation by 
minimization of the indicators I and J is a symmetric riskless full allocation. It satishes also the 
properties of comonotonic additivity, positive homogeneity, translation invariance, monotonicity, 
and continuity. 


We may write the indicators as follows: 

d 


and. 


/(«!, . . . , Ud) — ^ E ^ 

d , 

k=l ^ 

J (ui, . . . , Mrf) ^ ) E |ll|/jfc^Q|Il|^ii ^ 

d , 

~ X! ~ '^kY{Xk>Uk}^s^<>- 

k=l ^ 


Xi<u} 


E (^{Xk — Mfc)’^l{5<-u}) 5 


k=l 


X,>u} 


^ ® [i^k — Mfc)’^l{S>„} j . 

k=l 


In the special case where gk{x) = |a:|, O becomes: 

(1.2) P (Xj > Ui, S < u) = F {Xj > Uj, S < u) , V(i, j) G {1, 2,..., dY- 
For the J indicator, this condition can be written: 

(1.3) P (Xj > Wj, S' > u) = P {Xj > Uj, S > u) , V(i, j) G {1, 2,..., dY- 
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Some explicit and semi-explicit formulas for the optimal allocation can be obtained with this 
optimality condition. Our problem reduces to the study of this allocation depending on the nature 
of the distributions of the risk and on the form of dependence between them. 


2. Some general results in the independence case 

In this section we generalize the results presented in dimension 2 by Cenac et al. (2014) in the 
hrst section of their paper [S] to higher dimension. 

The results presented here give explicit forms to the optimal allocation for some specihc distribu¬ 
tions. This could be used as a benchmark to test optimization algorithms convergence. 


We also get some asymptotic results, when the capital u goes to inhnity. We study both the 
exponential and the sub-exponential cases, and we determine the difference between their asymp¬ 
totic behavior for exponential and Pareto distributions cases. 


We consider from now on that the penalty functions are identical and equal to the severity of 
local ruin gk{x) = g{x) = \x\,'ik G {1,2,... ,d}. The optimality conditions for minimizing the 
multivariate risk indicators I and J are given respectively by Equations (1.2) and (1.3). In this 
section we focus on the independence case. 


Recall that a* = ^ G [0,1], so that when u —)■ -f-cx), we may consider convergent subsequences 
in the proofs below. By abuse of notation, we consider lim at. In fact, we consider a convergent 

W^ + OO 

subsequence and get the existence of the limit by obtaining the uniqueness of the limit point. 


2.1. Independent exponentials. 

Assume Xi,X 2 ,... ,Xd are independent exponential random variables with respective parameters 
0 < Pi < /32 <■■■< /3d- Remark that in the particular case where Pi = P 2 = ■ ■ ■ = Pd, the optimal 
allocation is Ui = ■■■ = Ud = u/d. 


Proposition 2.1 (The optimal allocation for the indicator J). The allocation minimizing the risk 
indicator I is the unique solution in U^, of the following equations system: 

( 2 . 1 ) 

d 

h{Piai) - h{Pjaj) - ' (A - Z^^)) - h{o^j ■ {Pj - Pi))] = 0,V(i, j) G {1,2,.. .,dY, 

e=i 


where h is the function defined by h{x) = exp{—u-x), and denotes the constants Ai 
fori=l,...,d. 


d 


n 

j=i,j7^e 


A 

Pi-Pi’ 
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Proof. If Xi ~ £{Pi) are independent exponential random variables, then S' * = have a 

j=i\j¥=i 

generalized Erlang distribution with parameters {(3i, (32, ■ ■ ■, (3i-i, A+i) • • • > Pd), so we write: 


P I Xj > Mj, ^ Xj < M j = P (Xj > Mi) - P I Xj > Mi, ^ Xj > M 

J=1 


= FxPui) - FxPu) 


i=i 

Fs-Pu - s)fxAs)ds 


d d 

= h{(3iai) - h{(3i) - At>h{l3e)h{ai ■ (A - (Si)) + A(.h{(3i) 

i=i i=i 

d 

= h{/3iai) -J2Aih{Pe)h{ai ■ {Pi - Pi)), 


1=1 


because, Fx^x) = e Fs-i{x) = ^ 


d 

n 


p 


3 ^g-P^a; 


and '^Ai = l. 


i=l,li=-i j=l,j^f-,j^i ^3 Pi- 1=1 

The survival function of the generalized Erlang distribution with parameters {Pi, P 2 , ■ ■ ■, Pd) is 
given by: 

Pi 


cvw=i:( n . 

i=i j=i,j^i ^3 yi 


3 


= Y. 


1=1 


The optimal allocation is the unique solution in U!f, of the following equations system: 

P (xi > Mi, ^ Xfc < M^ = P (xj > Uj,Y^k < M ) , V(f, j) e {1,2,.. .,dy, 

\ k=l ) \ i=k ] 


which leads to (2.1). 


□ 


The resulting system is a system of nonlinear equations, which can be solved numerically. 

Proposition 2.2 (The asymptotic optimal allocation for the indicator /). When the capital u goes 
to infinity, the asymptotic optimal allocation satisfies: 

( 1 \ 


/Ml M2 Ud\ 

Im — = 

n^oo \u U M / 


Pi 


d 

.. 


Proof. Equations system (2.1) is equivalent to: 

(2.2) V(h j) e (1,2,... ,4^ h{Piai)[l-YAeh{{l - a^Pe)] = h{Pjaj)[l-Y Aeh{{l - afiPe)]. 

£=1 1=1 

Firstly, remark that for all i E (1,2,..., d}, lim sup ^ < 1 because if this result was not satished 
then there would exist i G {l,...,(i} such that lim— = lim Oi = 1, taking if necessary a 


U^OO 
U-j 


convergent subsequence of Oi. For all j ^ i, lim = lima,- = 0, and Equations system (1.2) 

U^OQ ^ U—^OO _ 

cannot be satisfied in this case. 

Equations system (2.2) is equivalent to: 


V(f,j) G (1,2, ...,4^ h{l3iai - PjOj) = 


l-j:tiAeh{{l-aj)Pe) 

l-j:tiAih{{l-ai)Pi)' 


6 












the right side of the last equations system tends to 1 when u tends to cx), therefore, we deduce 
that lim h{l5iai — Pjaj) = 1 and consequently: 


V(*, j) e {1,2,.. .,dy, lim a* = 


then, for all i G (1, 2,..., d}: 


lim CTj 

U^OO 


1 


A 


2^j=l 


1 • 

A 


A 

A 


lim ctj . 


□ 


Remark 2.3. Based on the above result, we can conclude that asymptotically: 

• since A < A have Oj > Oj, this means that we allocate more capital to the most risky 
business line. 

• ttj is a decreasing function of A- This observation is consistent with the previous conclusion. 

• aj is an increasing function of A for j A *• 


Proposition 2.4 (The optimal allocation for the indicator J). The allocation minimizing the risk 
indicator J is the unique solution in of the following equations system: 


(2.3) V(i,j) e {l,2,...,d}A ^klA(A)A(«* • (A - A)) - • (A - A))] = 0- 


e=i 


Proof. The proof is similar to that of Proposition 2.1 


□ 


Proposition 2.5 (The asymptotic optimal allocation for the indicator J). When the capital u 
goes to infinity, the optimal allocation minimizing the risk indicator J is the following: 

lim — = 1 and lim — = 0 V? G {2, 3,..., dj. 

U^OO 11 It—>CO 11 ^ 


Proof. Equations system (2.3) is equivalent to the following one: 

V(i, j) e (1,2,... ,d}A • h((i -Oi) • (A - A)) = IIA • K<^j ■ (A - A) + A - A)- 

i=i 1=1 


If lim ^ < 1, as M goes to +cx) in the equations of the previous system for i = 1 we get, taking if 
necessary a convergent subsequence, 

(2.4) Vj G (2,3, ...,d}, JimUkl^ • h(aj • (A - A) + A - A) = 4li. 

“ ^£=1 

The hrst terms of these equations can be decomposed into three parts as follows: 

ii™, H A ■ h{aj ■ (A - A) + A - A) = 4li ■ ■ (A - A)) 

i=l 

+ H • lim h{aj ■ (A - A) + A - A) 

i=2 

d 

+ H A ■ \in^h{aj ■ (A - A) + A - A)- 

e=j+i 
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For all, j > 1, '^Ai- lim h{aj ■ (fij — Pe) + I3i — Pi) = 0, because Pe — Pi >0 and Pj — Pe > 0 


1=2 


for I G {2,3, Moreover, ^ An - \mih{aj ■ [Pj — PP) + P^ — Pi) = 0, because for all 

£ G (j + 1, j + 2,..., d}, (Uj • {Pj — Pi) + Pi — Pi = {Pi — Pj){l — aj) + Pj — /Si > 0. So that 2.4 leads 
to Vj G (1,..., d} lim h{aj ■ {Pj — Pi)) = 1. We deduce that Vj G (2, 3,..., d}; aj o{^). 

This contradicts the necessary condition: lim ctj = 1- D 


2.2. Some distributions of the sub-exponential family. 

In most cases of risk distributions, we cannot give explicit or semi-explicit optimal allocation 
formulas, the difficulty comes from the lack of a simple form of the risks sum S and its joint 
distribution with each risk X^. In this section, we present asymptotic results (as u goes to inhnity) 
for the optimal allocation in the case of some distributions of the sub-exponential family. In this 
way, we generalize results of [H] to higher dimension. 

We recall the sub-exponential distributions family dehnition, consisting in distributions of positive 
support, with a distribution function that satishes: 


F*^{x) 

F{x) 




2 , 


where F*^ is the convolution of F. 

In Asmussen (2000) P, it is proven that the sub-exponential distributions satisfy also the following 
relation, for all d E N*: 

F*<^{x) , 

F{x) ’ 

where F*'^ is the d^^convolution of F. 

We shall use the following theorem proved in Cenac and al. 


Theorem 2.6 (Sub-exponential distributions |H])- Let X be a random variable with sub-exponential 
distribution Fx, Y a random variable with support independent ofX, and {u,v) G such 

that: 

• there exists 0 < < /s ;2 < 1 such that for u large enough, ni < ^ < ri 2 , 

Then, 

, P(X >v, X + Y>u) 

lim-=-= 1. 

Fx{u) 

2.2.1. The asymptotic behavior. 

Here we examine the asymptotic behavior of the optimal allocation by minimizing the indicators 
I and J in the cases of some sub-exponential distributions. 

In what follows, (ui,..., Ud) denotes the optimal allocation of u associated to the risk indicator J. 

Theorem 2.7. Let (Xi,X 2 ,... ,Xd) be continuous positive and independent random variables. 
Assume that V(i, j) G (1, 2,..., d}^: 

( 1 ) FxAx) ^^^-^Q{Fx,{x)), 

( 2 ) FxPs) ^^:^°^o{FxPt)),ift = o{s), 

^For {x,y) G M^+, we shall denote y = E){x) if there exist 0 < Ci < C 2 < oo, such that for x large enough, 
Ci<^<C2 

1 - J. - Z 












then, there exist Ki > 0 and K 2 < 1 such that, 


(2.5) 


Ki < — < ^2 G {1, 2,..., d}. 
u 


Note that the hrst condition of Theorem 2.7| is not satished for exponential distributions. How¬ 


ever, Pareto distributions satisfy the hypothesis of Theorem 2.7 


Proof. Taking if necessary a convergent subsequence, we assume that 3i G {1,... ,(i} such that: 
^ 1 or ^ implies that foil all j ^ i, ^ g^ then, it is sufficient 

to prove that the existence of an i G {1,..., d} such that — 0 is impossible. 


Let us assume the existence of i G {1,..., d}, such that: ^ g_ 

Then, 3j G {1,..., d}\i such that lim ^ g] 0, 1], therefore, Uj 
Using Assumptions (1) and (2), we deduce that: 


-|-cx) and Ui o{uj). 


( 2 . 6 ) 


M —> + 0 0 

Fx,{ui) ^ 


0 . 


The optimality condition (1.2) can also be written for all j ^ i as follows: 

(2.7) Fxi {ui) - P (Xj > Ui, S > u) = Fxj (uj) -F {Xi > Uj, S > u). 

That presents a trivial contradiction if Ui remains bounded. 

d 

Now, assume that Ui —)■ -t-cx). Recall that S~^ = ^ X^, then: 

k=l^k^i 


o{FxXu^)). 


P {Xi >Ui,S>u) = F (Xi > Ui, S-^ > s/u,S > u) +F (Xi > Ui, <s/u,S>u). 
We have: 

P (Xi > Ui, S-^ > ^/u,S > u) <F {Xi > Ui) P ( 5 -* > ^/u) 

Using assumption (2) and since Ui = o{u), 

P (Xj > Ui, S~^ < x/u. S' > m) < Fx,{u - y/u) 0{Fxi{Ui)) 

We deduce that: 

(2.8) 

We remark also that: 


P(Xi > Ui,S > u) o{Fxi{ui)). 


(2.9) P(X,- > u^,S> u) 0{Fx,{u,)) o{FxAui)). 


Equation (2.7) leads to: 


P (Xj > Ui,S > u) Fxj {uj) P {Xj > Uj, S > u) 


FxAui) 


Ti 


Fx_J^i) 

T 2 


FxAui) 


T 3 


Now, relations: (2.8), (2.6), and (2.9), imply that Ti, T 2 , and T 3 , go to zero, and this is a contra¬ 
diction. □ 

Proposition 2.8. Let Xi,... ,Xrf be continuous, positive and independent random variables such 
that the support of the density of (Xj, S) is (M+)^. Let {ui,... ,Ud) be the optimal allocation of u 
associated to the risk indicator F We assume: 
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( 1 ) there exist 0 < Ki < K 2 < 1 such that for all i = 1,..., d and for all u G 


Ki < < K2, 

U 


( 2 ) for all i = 1,... ,d, if y = y{x) is such that 

0 < lim inf — < lim sup — < 1 

x—¥oo nr 


then 


■^<limsup^ 

X X —>00 X 

x^op 


Then, for all i,j = l,...,d, 


FxAy) 

FxA'^i) 


0 . 


1 . 


F Xj{uj) 

Assumptions of Proposition |2.8| are satisfied for distributions of exponential type (see remark (2.9) 
below). Its application gives another proof to Proposition 2.2 In contrast, Proposition 2.8 cannot 
be used for Pareto distributions. 

Proof. Following the lines of the proof of Theorem 2T, take 0 < 7 < 1 — ^ 2 , 

P(Xj > Ui, S > u) = P(Xj > Ui, S~^ > 7M, S > u) + P(Xj > Ui, S~^ > 'ju, S > u). 

As before, 

P(Xi > Ui, S~" > lu, S > u) < P(Xi > Mi)P(S'“* > 7 m) o{Fxi{ui)). 

On the other hand, 

P(Xj > Ui, S~^ > 'ju, S > u) < P(Xj > u — 7m) 

Ki 'll • 

because 0 < -< 


1 — 7 (1 — 7 )m 1 


<^< 1 . 


■7 


So that, P(Xj > Ui, S > u) o{FXi{ui)) and the same computation gives P(Xj > Uj, S > 

u) o{Fxj{uj)). Now, Ui and Uj satisfy Equation (1) and thus. 


/ X >-i-oo Fx-iuj) , ^ FX (uj) 

1 +0(1) =777+o 1)77=7. 

FxA'Ui) FxAti,) 

^ X iX' ) 

This implies that \ \ is bounded from above and thus 

1 + 0 ( 1 ). 

FxAUi) 


□ 


Remark 2.9. We remark that the hypothesis of Proposition 2^ are satished for distribution of 
exponential type, that is distributions verifying: 


F^^(x) = 0 (e-^'^). 
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Indeed, in this case, we have 0 < lim — < lim — < 1, Vi G {1,..., d}. In fact, if Ui o(m). 


then, 3j G {1,... such that ^ t k G]0, 1], so Ui o{uj). 


Since > 0, i^tUi o{fijUj). As in the proof of Theorem 2.7, we get: 
and. 

From the optimality condition. 


P {Xi >Ui,S> u) o{FxM)) 


P {Xj > Uj, S > u) o{Fxj{uj)). 


Fxi{ui) P (Xj > Ui,S > u) _ Fx^{uj) P {Xj > Uj, S > u) 








^3 


which is absurd because as m —)■ +cxo: 

. FxM) = 

• Ti = o(l), 

• Ta = sin(;;e = oifijUj), 

• and T 3 = —)• 0 . 

Proposition 2.10 (The asymptotic optimal allocation for the indicator I). Under the same con¬ 
ditions of Theorem 2.7, and if, for all i G {!,... ,d}, Fxi is a sub-exponential distribution, that 
verifies: 

Fj^ 0 ( 1 ), for 0 < < - < K 2 < 1 . 

Fx,{x) X 

Then, by minimizing the I indicator, Ui and Uj satisfy: 

( 2 . 10 ) FxM) - Fx,{u) FxM,) - FxM) + o{FxAn)). 


Proposition 2.10 is applicable in the case of Pareto distributions. So, we will use it for deter¬ 
mining the optimal asymptotic allocation, for independent risks of Pareto distributions in the next 
subsection. 


Proof. The proof of this theorem is a direct application of Theorems 2.6 and 2.7 


□ 


Now, we focus on the asymptotic optimal allocation by minimizing the risk indicator J, and we 
study the case of sub-exponential distribution family. 

Proposition 2.11 (The asymptotic optimal allocation for the indicator J). Let (Xi,X 2 ,... ,Xd) 
be continuous positive and independent random variables, such that there exists i G {l,...,(i} 
with a sub-exponential distribution, the optimal capital allocation by minimizing the J indicator 
(til,... ,Ud) verifies, for all j 7 ^ i: 

, P(X,' > Uj, S > u) 

lim —^^^ = 1 . 

Fxiiu) 


Proof. The solution to (|1.3) satishes : 


VjG{l,2,...,d}, 


P {Xj > Uj, S > u) P {Xj > Uj,S > u)) 


P {Xj >u) P {Xj > u) 

When u goes to -fcxo, and using Theorem |2.6[ we obtain Proposition 2.11[ 


□ 
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2 . 2 . 2 . Application to Pareto independent distributions. 

We consider d independent random variables (Xi, X 2 ,..., X^) of Pareto distributions, with param¬ 
eters (a, &i){i=i, 2 ,...,d} respectively, such that bi > b 2 > ■ ■ ■ > bd > 0. Therefore, these distributions 
will be characterized by densities and survival functions of the following forms: 

— a—1 




and 


Fx.W = (l + r 


X 


Proposition 2.12 (The asymptotic optimal allocation minimizing I). Asymptotically, the unique 
solution to (1.2) satisfies: 


V(i, j) e {1,2,.. .,dy 


( lim Oj 

I u—^co 


( lim ttj 
I U—^CXD 


Proof. Follows from Proposition 2.10 


□ 


Proposition 2.13 (The asymptotic optimal allocation minimizing J). The unique solution to (1.3) 
satisfies: 

lim oi = 1 and lim Oj = 0, Vi G |2, 3,..., dj. 


Proof. We suppose that 3 G (1,..., d} such that 0 < lim aj < 1. 

From Theorem 12.61 : 

, IP(X,' > Uj, S > u) 
lim - = 1 . 

Fxfiu) 

On the other hand, and applying Proposition 2.11| in the Pareto distributions case, we get for 
i G (1,... ,d} \ {j}: 

P(Xj >Uj,S>u) _ P(Xj >Uj,S>u) Fxfiu) u^+oc Fxfiu) _ A + ^ 


Fx^u) 

then, for i G (1,..., d} \ {j}: 


Fxfiu) 


FxAu) FxAu) 


l + F 


p(x,>^„^>^) ^ 

u^oo yjj. J 

That is absurd. We deduce that Vi G (1,2,..., d}: 

lim Qfj G (0,1}, 


and since I]f=i cti = 1, then, there is a unique i such that lim = 1 and for all j A i lim aj = 0 . 
Let us recall the dehnition of the order stochastic dominance, as it is presented in Shaked and 
Shanthikumar (2007) [SO]- For random variables X and Y, Y hrst-order stochastically dominates 
X if and only if: 

Fx{x) < Fy{x), 'ix G MA, 
and in this case we denote: X =4st Y. 

Now, clearly X^ ^st ■ ■ ■ Ast X 2 ^st Xi because bi > b 2 > ■ ■ ■ > bd > 0. 

In [T7] (Proposition 3.10) we have proved that the optimal capital allocation satishes the mono¬ 
tonicity property. We deduce that ai > ■ ■ ■ > ad, and thus lim ai = 1. □ 

11^00 
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2.3. Comparison between the asymptotic behaviors in exponential and sub-exponential 
cases. Let us consider the following risk indicator: 

d d 

hoc{Ui^ . . . , Ud) = ^ E ({Xk — = X! ((^fc ~ . 

k=l k=l 


With this indicator, the impact of each branch on the group is not taken into account. This 
indicator thus takes only local effects into account. If, for any k = 1,... ,d, the random vector 
(Xfc, S') admits a density whose support contains [0 ,m]^, then the optimality condition associated 
to hoc gives: 

P(Xj > Ui) = P(Wj > Uj). 


We remark that this corresponds to the asymptotic result for indicator I of Proposition |2.8[ where 

^ other words, in the case of exponential type and independent risks, the 

group effect is asymptotically negligible. This behavior is also clear in Proposition 2.2[ where we 
found the asymptotic capital allocation for I for independent exponential distributions as: 


Ui = — . u, for alH = 1 ,..., d, 
-L 


_ 


as for hoc- 


Proposition 2.10 shows that for indepe ndent sub-exponential distributions the asymptotic be¬ 
havior is different. Indeed, in Equation (2T0|, the terms Fxi{u) and Fxj{u) lead to take into 
account the group effect. In the Pareto distributions case as example, recall that Proposition 
gives that the asymptotic behavior of the optimal allocation for I is described by : 


2.12 


lim Uj 

U^OO 


lim a, 

U^OO 


whereas the optimal allocation for hoc for independent Pareto distributions is given by ai 
for alH G {1 ,..., d}. 


bj 


For optimal allocation by minimizing the J indicator, the asymptotic behavior is identical for 
the two families of distributions, the riskiest branch is considered as hrst responsible of the overall 
ruin, and thus, the optimal solution is to allocate the entire capital u to this business line. 


3. The impact of the dependence structure 

In this section, we focus on the impact of the dependence structure on the optimal allocation. 
We study at hrst the impact of mixture exponential-gamma to construct a correlated Pareto 
distributions, compared to the independence case presented in the previous section. Then we 
analyze the optimal allocation in the case of comonotonic risks. The last sub-section is devoted to 
the study of the impact of the dependence nature on the optimal allocation, using some bivariate 
models with copulas. 
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3.1. Correlated Pareto. 

Let {Xi,..., Xd) be a mixture of exponential distributions such that for all i G {1,2,... ,d}, 
Xi ~ S{Pi9), with (/5i < /52 • • ■ < Pd), and 9 ~ r(a, 6 ). Therefore, Xi have survival functions of the 
form: 

FxAx) = Fx^\e=efe{9)d9 = fe{9)d9 = (^1 + , 

consequently, X^ have Pareto distribution of parameters (^a, They are conditionally indepen¬ 
dents. So, the idea is conditioning on the random variable 9 and then integrate the formulas found 
for the case of independent exponential distributions. This model has been studied in e.g. mn 

Recall Ai= n , i = !,■■■, d. 

j=i,j^e Pi 

Proposition 3.1 (The optimal allocation for the indicator J). The optimal allocation minimizing 
the multivariate risk indicator I is the unique solution in of the following equations system: 


(3.1) 


s{/3iai) - si/SjOj 


d 

E 

i=i 


Ae[s{aiPi -|- (1 — apPi)) — s{ajPj + (1 


^j)Pe))] = 0 , 


where s is the function defined by s{x) = (1 -|- x|) “ and OLi = — for all i E [I,... ,d}. 


Proof. It suffices to integrate Equations System (2.1), multiplied by the density function of 9. □ 


Proposition 3.2 (The asymptotic optimal allocation for the indicator /). When the capital u goes 
to infinity, the optimal allocation by minimization of the risk indicator I is the unique solution in 
of the following equations system: 


W{i,j) e {l,2,...,dy. 


(3.2) 


(Piai) “ - iPjaf “ - E Mi<^iPi + (1 - cti)Pe) “ - {ajPj + (1 - afPe) “] = 0. 


e=i 


Proof. We divide Equations System (3.1) by s(l), and let u go to -|-cx) to get Equations System (3.2). 

□ 


Proposition 3.2 shows the impact of the dependence related to the mixture. Indeed, in the case 
of independent Pareto distributions, of parameters (a, J-) , the asymptotic optimal allocation 

for the indicator I is given by Proposition 2.12 as the solution of the equations system: 

V(q j) E [1,2,..., = (A)-“ - (/d,)-“ . 

Each equation in this system depends only on two risks, unlike the mixture case, where the 
equations of Equations System (3.2), depend on all the risks. 

Proposition 3.3 (The optimal allocation for the indicator J). The optimal allocation minimizing 
the multivariate risk indicator J is the unique solution in of the following equations system: 


(3.3) V(i, j) e {1,2,.. .,dy,J2Msi(^iPi + (1 - O'PPe)) - s{ajPj + (1 - afPt))] = 0. 

i=i 


Proof. It suffices to integrate Equations System (2.3), multiplied by the density function of 9. □ 
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Proposition 3.4 (The asymptotic optimal allocation for the indicator J). When the capital u goes 
to infinity, the optimal allocation by minimization of the risk indicator J, is the unique solution in 
of the following equations system: 

d 

(3.4) V(i, j) G {1,2,.. .,dy, + (1 - + (1 “ = 0. 

£=1 


Proof. We divide Equations system (3.3) by s(l), and we let u fo to +cxd to get Equations sys¬ 
tem (3.4). □ 


Proposition 3^ shows that for the indicator J, the asymptotic behaviour of the optimal capital 
allocation takes into account the mixture effect. In fact, for independent Pareto distributions, we 


have proved in Proposition 2.13 that we allocate the entire capital u to the riskiest branch Xi, 
while the asymptotic optimal allocation in the correlated Pareto distributions case is the solution 


of Equations system (3.4). 


3.2. Comonotonic risks. 

The concept of comonotonic random variables is related to the studies of Hoeffding (1940) [13] 
and Frechet (1951) [12]. Here we use the dehnition of comonotonic risks as it was hrst mentioned 
in the actuarial literature in Borch (1962) [3]. 

A vector of random variables (Xi, X 2 ,..., X„) is comonotonic if and only if there exists a random 
variable Y and non-decreasing functions (pi,..., (p„ such that: 

(Xi,...,x„)^((pi(x),...,(p„(y)). 

In the case where the risks Xi,... ,Xd are comonotonic, we may give explicit formulas for the 
optimal allocation minimizing the multivariate risk indicators I and J, and for some risk models. 
For that, we use the existence of a uniform random variable U such that Xj = Fx^iU) for all 
i e (1,... ,d}, and S = Yfi=i Fx^iU) = (p{U), where ip{t) = Ef=i Fx] (t), 93 is a non-decreasing 
function. 

The main result of this section is given below. 


Proposition 3.5. Let Xi,... ,X^ he comonotonic risks, with increasing distribution functions and 
support containing [0,m]. The optimal allocations for indicators I, J and fioc coincide, they are 
given by {ui,, Ud) G and 

Fxfiui) = Fxfiuj) y i,j = l,...,d. 

Proof. Let us denote: Wi = Fxfiui), v = ip~^{u), Mi = max(t(;j,u). The indicators / and J may 
be rewritten for {u\,..., Ud) G U!f\ 

d 

/(wi, . . . , tic/) ^ ^ IE (^) U 

i=l 
d 

.J{Ui,...,Ud) = -Ui)l{u>Mi}) ■ 

i=l 

We remark that since (ui,... ,Ud) G and Fx^ is strictly increasing for alH G (1,..., d}, we 
cannot have that Wi < v for all i, so that, I is not trivially equal to 0. We use Lagrange multiplier 
to get that the minimum of I and J are reached in respectively for: 

• F{U > Wi, U <v) = F{U > Wj, U < v), for i, j = 1,... ,d, 

• F{U > Mi) = F{U > Mj), for i,j = l,...,d. 

These equality are acheaved if and only if Wi = Wj = v or in other words if Fx.iu.) = F,cAuj). We 
remark that the minimum of I is then 0. □ 
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The following three corollaries are direct applications of Proposition 3.5 to some particular cases. 


Example 3.6 (Comonotonic exponential model). For comonotonic risks of exponential distribu¬ 
tions Xi ~ exp(/3j), the optimal allocation by minimization of the two risk indicators is: 


Vi e {1,.. .,d}, Ui 


1/A 

—^— u- 

eU va 


It is noticeable that in this particular case, the optimal allocation for comonotonic risks coincide 
with the asymptotic allocation in the independence case (see Proposition |2.2[). 


Example 3.7 (Comonotonic log-normal model). For comonotonic risks of log-normal distributions 
Xi ~ LN{fii,a), the optimal allocation by minimization of the two risk indicators is: 


Vi G {1,... ,d}, Ui 


exp(pi) 

, , , u. 

J2e=i exp(p£) 


Example 3.8 (Comonotonic Pareto model). For comonotonic risks of Pareto distributions of the 
same shape parameter a: Xi ~ Pa{a,Xi), the optimal allocation by minimization of the two risk 
indicators is: 

Vi e {1,... ,d}, Ui = , \ u. 

Xi=i M 


Contrary to the exponential case, this result does not coincide with the asymptotic behavior 
obtained for independent Pareto distributions (see Proposition 2.12). 


3.3. The dependence impact with some copulas models. 

In this section, we study the impact of the dependence on the optimal capital allocation using 
some copulas (see Nelsen iia for review on copulas). The idea is to hnd the optimal allocation as 
a function of the copula parameters in each case. We focus on the indicator I, the same kind of 
calculations can be done for the indicator J. 


3.3.1. FGM Bivariate Model. 

Let Xi and X 2 be two risks of marginal exponential distributions Xi ~ exp(/5j) and FGM bivariate 
dependence structure with —1 < 0 < 1 as parameter (see Nelsen [T5] .Example 3.12., section 3.2.5). 
We assume that A < A/2- 

In this case, the copula Pearson correlation coefficient is given by pp = |, and the bivariate 
distribution function is: 

Fx„x,{xi,X2) = (1 - e-F^^){l - + ^(1 " 


Proposition 3.9 (The optimal capital allocation for the indicator I in the FGM Model). For the 
indicator I, the optimal allocation of a capital u is given by (/3u, (1 — /3)u) such that (3 = ui/u is 
the unigue solution in [0,1] of the following eguation: 

(3.5) (1 + 2A(h(A - h{a - a/3)) + 2e{h{2/3) - h{2a - 2a/3)) 

-l-(l -I- 9)h{a 3- /3 — a/3) -|- 6h{2a 3-2/3 — 2a/3)—9h{a 3-2/3 — a/3) — 9h{2a 3- /3 — 2a/3) 

= — ^^{h{a) 3- ah{l)) 3 - {h(2a) 3- ah{2)) - - — {2h{a) -3- ah{2)) - {h{2a) -3- 2ah{l)), 

a — 1 a — 1 a — 2 2 a — 1 

where, h is the function h{x) = exp{—/3iux),and a = A/A- 


Proof. The proof is postponed to Appendix A.l 


□ 


Remark 3.10. In the case of 0 = 0, we find exactly Equation (2.1) given by Proposition 2.1 for the 
independent exponential distributions model. 
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Equation (3.5) gives the behavior of the optimal allocation with respect to 9. It may be solved 
numerically. 

Figure [^presents an illustration of the optimal allocation variation with respect to the dependence 
parameter of the FGM copula. 


The optimal allocation according to the dependence parameter In the FGM model 



The dependence parameter 0 


Figure 1. /? as a function of 9. Case : /^i = 0.05, (32 = 0.25 , and u = 50 

For the illustration parameters, we remark that (3 is an increasing function of 9, this can be 
verihed analytically using the implicit function theorem. It is important to remark that the optimal 
allocation is also a function of the capital u. The variations of (3 with respect to the dependence 
parameter depends on u. Figure]^ give what we obtain as result for an allocation with the same 
distributions parameters but with u = 100 as allocation capital. 


The optimal allocation according to the dependence parameter in the FGM model 



Figure 2. /3 as a function of 9. Case : j3i = 0.05, (32 = 0.25 , and u = 100 

The optimal allocation depends on the initial capital. As an example, with the same parameters, 
/5 is a decreasing function of 6* if u < 53, and it is an increasing function of 9 for u > 53. 

The variation range size of [3 is very small, and its size is function of the distributions parameters 
and the allocation capital. That is due to the dependence structure, the FCM copula present only 
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a light dependence. For Clayton copula as example, which is stronger dependence structure, the 
variation range size of (3 is more important. 


3.3.2. Marshall-Olkin Model. 

Let Yi ~ exp(Aj), with z = 0,1, 2 be three independent random variables. 

We construct two random variables with common shock: Xj = min(yj, Yq) for z = 1, 2. Xj’s have 
exponential marginal distributions of parameters /S* = A* + Aq (see e.g. Nelsen |TH] section 3.1.1.). 
This dependence construction model has as Pearson correlation coefficient pp = 

The joint distribution function is given by: 


Fxi,X 2 {xi,X 2 ) = P(Xi > Xi,X 2 > X 2 ) = P(Ti > xi,Y 2 > X 2 ,Yo > max(a;i,a; 2 )) 

_ g-Aia;ig-A2a;2g-Ao max(a:i,a;2) 

_ g-(Ao+Ai)a:ig-(Ao+A2)a:2gAo min(a:i,X2) 


and the joint density function is the following: 


( fxi,X 2 (^hX 2 ) = /3ie F^^(/32-Ao)e si xi>X 2 

fx^,X2{xi,X2) = < fxi,X2i^l^^2) = {(3l - Xo)e-^hi-^o)xip^^-(^2X2 gi Xi <X2 

{ /xi,X 2 (^i’^2) = Aoe^^i'^e^^^a^gAo* gj 2;;^ = X2 = a; 

Proposition 3.11 (The optimal capital allocation for the indicator I in the Marshall-Olkin 

Model). fTe suppose that Ai < A 2 . The optimal allocation of a capital u minimizing the indicator 
I is given by {(3u, (1 — (3)u), such that {3 = ui/u is the unigue solution in [0,1] of the following 
eguation: 


5 '(/^ 2 (l — 13)) — g{l3il3) + g{{/3i — ^ 2)13 + A 2 ) + - 5 —^- 5'((-^2 — /5i)(l — /3) + /3i) 


^ ' /3,-X2 

-9(/?2) = XV9(A)+9(V2)1 


A, 


Ai —132' 


(3i — X 2 


I3i — X 2 Ai — /32 ’ 


where, A^ = Aq + Ai + A 2 , and g is the function g{x) = exp{—ux). 


Proof. The proof is postponed to Appendix |A.2[ 


□ 


Remark 3.12. In the case Aq = 0 which is the independence case, we hnd again Equation (2.1) 


given by Proposition 2.1 


We can consider Aq as a dependence parameter in this model. Figure [^presents an illustration 
of the variation of /3 as a function of Aq. 
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The optimal allocation according to the dependence parameter of the MO Model 



The dependence parameter Xo 


Figure 3. /3 as a function of Aq : Pi = 0.05, P 2 = 0.25 and u = 50. 

One can notice that /? is a decreasing function of Aq. This is coherent with the increase of P as 
a function of a = ^ = demonstrated in [S] in the independence case, since the two risks are 

independent, conditionally to Yq. 


Conclusion 

In this article, we have studied the allocation asymptotic behavior based on the level of the 
group capital. It has enabled us to build an idea of the capital level impact on the sensitivity of its 
allocation between branches. The comparison between the asymptotic optimal allocation, in the 
case of sub-exponential and exponential distributions, underscores the impact of the risks nature 
on the behavior of the allocation for a very large capital. 

Compared to the classical risk allocation methods, the allocation by minimizing multivariate 
risk indicator take into account the capital level. This seems to be a more acceptable behavior in 
a capital allocation operation. 

The risk aggregation transforms the global risk portfolio in a univariate risk. At the opposite, 
the allocation is based on a multivariate analysis. The goal of this paper is to highlight the 
importance of the dependence modeling in the success of an allocation capital operation. The 
capital allocation is also sensitive to the risks nature, we tried in this paper to make more evident 
the impact of distributions’ nature on the allocation behavior. Since the capital allocation is an 
important hnancial decision, especially for groups, insurers must be very careful in their risk and 
dependence modeling choices to get an efficient allocation. 
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Appendix A. Proofs 


A.l. Proof of Proposition [37^ 

Proof. The bivariate density function is the following: 

fxi,X 2 {Xl,X 2 ) = (1 + 0)f{Xi,X2,/3i,/32) + ef{xi,X2,2/3i,2/32) 

- 9f{xi,X2,2(3i,(32) - 9f {xi,X 2 , (3i, 2(32) fxr,X 2 {xi,X 2 ), 
where / is the function f{x, t, a, b) = 

We use the equality: fxi,s=Xi+X 2 {xi, s) = fxi,X 2 {xi,s - a;i)l{s>xi} for all s > xi, to hnd the 
expression of Fx,, s(a;i,s), using a double integration: 

PS PX\ PX\ PS 

Fxi,sixi,s)= / fxi,X 2 ix,t - x)l{t>x}dxdt = / fxi,X 2 {x,t-x)dtdx 

Jo Jo Jo Jx 

= (1 J- 9)F{xi, s,l3i, (32) + 9F[xi, s, 2(3i, 2(32) — 9F{xi, s, 2(31,(32) — 9F{xi, s, (3i, 2(32), 
where F is the following function: 

F(xu s, a, b) = r r abeF^-^^^e-^^dtdx = 1 - - ^_^-bs+(b-a)x, _ 

Jo Jx b — a b — a 

The same way and by the symmetry of the FGM model: 

Fx 2 ,six 2 , s) = (1 -b 9 )F{x2 , s, (32, (3i) + 9F{x2, s, 2(32,2(3i) - 9F{x2, s, 2(32, (3i) - 9F{x2, s, (32,2(3i). 

Using P(Xj > Ui, S < u) = P(S' < u) — P(Xj < Ui,S < u), the optimal allocation is the 
unique solution in of the equation: Fxi,s(wi,m) = Fx 2 ,siu 2 ,u). Then, the optimal allocation is 
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determined by (3 the solution of the equation: = Fx 2 ,s{i^ ~ (3)u,u). 

Since, 

Fx,,s{Pu, u) = l + 49-{l + 29)h{P) - 29h{2^) 

+ (1 + 9) - [h{a) — h{a + /3 — a/3)] + 9 -[h(2Q;) — h{2a + 2/3 — 2a/3)] 

a — 1 a — 1 

— 9 - [h{a) — h{a + 2/3 — a/3)] — 9 - [h{2a) — h(2a + (3 — 2a/3)], 

a — 2 2a — 1 

and, 


i^X 2 ,s((l - I3)u, u) = 1 + 40 - (1 + 20)h(a(l - /3)) - 29h{2a{l - /3)) 

CY (Y 

+ (1 + 9)- -[h(l) — h{a + (3 — a/?)] + 9- -[h(2) — h{2a + 2/3 — 2a/3)] 

1 — a 1 — a 

^(Y (Y 

— 9- - —[h(l) — h{2a + /3 — 2a/3)] — 9- - [h{2) — h{a + 2/3 — ct/S)], 

1 — 2a 2 — a 


we deduce from that the equation presented in the proposition 3.9 


□ 


A. 2. Proof of Proposition 3.11 


Proof. The joint distribution function is given by: 

rs pXi 


rs px\ 

Fx,,s{xi,s)= / fxus{x,t - x)l{t>x}dxdt 
Jo Jo 

PS PX\ PS PX\ 

= ^ ~ x)l{2x>t>x}dxdt + ^ ~ x)l{t>2x}dxdt 

PS PXl 

+ - ^)^{t=2x}dxdt. 


we distinguish between two cases: 
Case s > 2xi: in this case. 


PS PXl PS PXl 

^ - x)l{2x>t>x}dxdt = fxi,si^y t - x)l{2x>t>x}dxdt 

p2xi PXl 


dxdt 


1 -2x1 r^i 

fx^^si^yt - x)l{2x>t>x}dxdt 


/o ^0 

p7nin(xi,t) 

Jo Jt /2 

PXl pt 


fxi ^ “ x)dxdt 

2xi PXl 

xi Jt/2 


PXl Pt pZXi PXl 

/ / fxusi^^t-^)dxdt+ / / fx.si.Xyt - x)dxdt, 

Jo Jt/2 Jx, Jt/2 


and, 

f'S PXl 


PS PXl PXl PS PXl PS 

fx^,si^^'t - x)l{t>2x}dxdt = fx^ six^t - x)l{t>2x}dtdx = fx^ six^t - x)dtdx, 

and. 


PS PXl 

fxi,si^^ t - x)l{t=2x}dxdt = 


Ar 


Aq + Ai + A2 


_ p— (Ao+Ai+A2)3^1 


(1-e 


), 
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then, we deduce the explicit expression of Fx, 


Fx^,s{xi,s) = 


2/3iA2 


(/9i-A2)(/3i + A2 ) 
( 1 -e 




( 1 -e 


Ao 




Pi — A 2 


/Si 


Pi — A2 


(g-/3ia:i _ g-(/3i+A2)xi^ 


'^1 _ ^-(Ai+/32)a:i'| _ 


+ 


Ai + P 2 

Ao 


Aq + Ai + A2 
Case 2x\ > s > xi ■. 

rs rxi 


Ai — P 2 

_ g —(Ao + Ai+A2)3^1^ 

s rmin{xi,t) 


I ^-(Ai-fe)a:i-/l 2 S 

Ai - /?2 


^5 prmn[xi,t) 

/xi.5(^) ^ - x)l{2x>t>x}dxdt = t - x)dxdt 


t/2 

X\ ft 


and, 


and. 


rxi rt rs rxi 

= / / fxusi^^t - ^)dxdt + / / fx,s{x,t - x)dxdt, 

JO Jt/2 ’ Jx^Jtl2 

fs fxi fs ft/2 

J J fxi,si^J-^)^{t> 2 x}dxdt = J J fx^ six,t - x)dxdt, 


f f ^ - x)l{t=2x}dxdt = , (1 - e Co+Ai+A2)s/2^^ 

Jo Jo An + Ai + Ao 


Aq + Ai + A 2 

then, we deduce also in this case, the explicit expression of Fxj^,s{xi, s): 


Fx^,s{xi,s) = 


2 P 1 X 2 


(/Si - A 2 )(/Si + A2) 




Ao 


Pi — A2 


(1 _g-/3ixi) 


(Si r ^ — ^-lX^ _ —f/li — Ao'lii — A2S^ I "^1 /-I 3—/32 s\ 


Pi — A2 

2Ai/32 


(g-^ia:i _ g-(/3i-A2)a;i-A2S^ _|_ 


Al — P 2 


(1 - 


-(1-e 


-(Ai+/32)s, 




An 


-( 1 -e 


~(Ao + Ai+A2)5/2\ 


{Xi — P2){Xi + P 2 ) ' Ao + Al + A 2 

We remark that Aq + Ai + A2 = Xi + P2 = A2 + Pi, and we suppose that Ai > A2. Using the 
monotony property, we deduce that 1 > /3 > 1 / 2 , then 2Pu > u > Pu. So, for ui = Pu, and 
g{x) = exp(— xm): 

Fx^,s{Pu,u) = l-g{PiP)--^J—g{P2) + -x^F^gi{Pi-X2)P + X2)+g{\s/2)[ ' 


Al — P2‘ 


Pi — Ao 


Al — P2 Pi — X2 ' 


and. 


Fx2,s{{^ - P)u, u) = l- g{P2{l - P) + 


Ag 


A 2 —Pi 


[ 5'((-^2 — Pi) 0 - ~ P) F Pi) ~ 9 {Pi)]- 


That is sufficient to get |3.11 


□ 
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